Effective pairing interaction in the two-dimensional Hubbard model 
within a spin rotationally invariant approach 
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We implement the rotationally-invariant formulation of the two-dimensional Hubbard model, 
with nearest-neighbors hopping t, which allows for the analytical study of the system in the low- 
energy limit. Both U(l) and SU(2) gauge transformations are used to factorize the charge and 
spin contributions to the original electron operator in terms of the corresponding gauge fields. The 
Hubbard Coulomb energy U term is then expressed in terms of quantum phase variables conjugate 
to the local charge and variable spin-quantization axis, providing a useful representation of strongly 
correlated systems. It is shown that these gauge fields play a similar role as phonons in the BCS 
theory: they act as the "glue" for fermion pairing. By tracing out gauge degrees of freedom, the 
form of paired states is established and the strength of the pairing potential is determined. It is 
found that the attractive pairing potential in the effective low-energy fermionic action is non-zero 
in a rather narrow range of U/t. 



I. INTRODUCTION 

Superconductivity represents a remarkable phe- 
nomenon where quantum coherence effects appear at 
macroscopic scalcii The quantum-mechanical phase of 
the electrons gains rigidity and as a result the prop- 
erties of the quantum wave show up at the macro- 
macroscopic level. Thus, the superconducting properties 
are the manifestations of the spontaneous breakdown of 
one of the fundamental symmetries of matter, namely, 
the U(l) gauge symmetry. The discovery of the cuprate 
superconductorsS has sparked a widespread interest in 
physics which goes beyond the traditional Fermi-liquid 
framework usually employed for understanding the ef- 
fect of interactions in metals. The question of whether 
the pairing interaction in the cuprate superconductors 
should be characterized as arising from a pairing glue has 
recently been raised.'^ While there is a growing consen- 
sus that superconductivity in the high- Tc cuprates arises 
from strong short-range Coulomb interactions between 
electrons rather than the traditional electron-phonon in- 
teraction, the precise nature of the pairing interaction re- 
mains controversial. In this context the Hubbard model 
is considered as essential physical system for treating su- 
perconductivity in the strongly correlated electron sys- 
tems and has been intensively studied with a variety of 
methods such as quantum Monte-Carlo^ (QMC), exact 
diagonalization,^'^ path-integral renormalization-group,® 
functional renormalization-group)^ and various quantum 
cluster methodsiiS As a principal model describing the 
electronic correlation in the system, the Hubbard model 
has been used in many works to study the pairing in- 
stabilities which as usual are given by the second-order 
effective interaction with respect to the Coulomb inter- 
action. In this context the structure of the pairing in- 
teraction, the two-dimensional (2D) Hubbard model, has 
been recently analyzed in— li^iH^ where the dynamical 



cluster Monte Carlo approximation is applied to two- 
dimensional Hubbard model with nearest-neighbors hop- 
ping and on site Coulomb interaction. The Monte Carlo 
simulations have been also employed to study the phase 
separation and pairing in the doped two-dimensional 
Hubbard model." However, the question whether the 
Hubbard model even supports superconductivity without 
additional interactions remains a subject of controversy. 
Different mean-field theories suggest conflicting ground- 
state order parameters and correlations, while finite-size 
QMC simulations for the doped 2D Hubbard model in 
the intermediate coupling regime of correlation energy U 
support in general the idea of a spin-fluctuation-driven 
interaction mediating d-wave superconductivity. How- 
ever, the fermion sign problem, limits these calculations 
to temperatures too high to study a possible transition. 
These simulations are also restricted to relatively small 
system sizes so that the off-diagonal long-range order 
has not been ascertained. For theoretical understand- 
ing of the mechanism of superconductivity in cuprates, 
the knowledge of bosons mediating the pairing is of piv- 
otal importance. Here, the underlying attraction force 
appears very puzzling since it is hard to reconcile the mi- 
croscopic attractive interaction with the completely re- 
pulsive bare electron-electron forces. This issue is closely 
related to the construction of the low-energy effective 
theory for the electronic system. A powerful tool for 
the quantitative investigation of microscopic models is 
provided by the study of effective theories: if one is 
able to single out the most relevant low-energy config- 
urations, an effective theory can be extracted from the 
microscopic lattice Hamiltonian. This procedure is often 
implemented via the projective transformation, which re- 
sults in removing of high-energy degrees of freedom and 
replacing them with kinematical constraints as exempli- 
fied, e.g., by the t-J modeli^ In such approaches, the 
high-energy scale associated with the charge gap is ar- 
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gued to be irrelevant, hence the focus exclusively on the 
spin sector to characterize the Mott insulator. However, 
the charge-transfer nature of the cuprates plays an es- 
sential role in the doped systems, so that with discard- 
ing charge degrees of freedom an important part of the 
physics may be lost. In the same spirit a detour from 
the strict projection program was recently proposed in 
a form of the "gossamer" superconductor,^^ recognizing 
the role of the double-occupancy charge configurations. 
However, the most interesting and relevant situation of 
strongly correlated systems, where magnetic as well as 
charge degrees of freedom interact, was until now inves- 
tigated to a much lesser extent since it requires the treat- 
ment of the Hubbard Hamiltonian without imposing any 
restrictions on the correlation energy U. 

In the present paper we study the emergence of the at- 
tractive pairing interaction in the two-dimensional Hub- 
bard model by resorting to the analytical method that is 
deeply rooted in the inherent spin-rotational and gauge- 
charge symmetries of the model. To keep the spin- 
rotationally invariance, we write the action of the sys- 
tem using other bosonic and fermionic variables which 
are introduced with appropriate U(l) and SU(2) trans- 
formations. We construct a SU(2) spin-rotational and 
charge U(l) invariant theory using the electron operator 
factorization. ^^'^^ Furthermore, we derive the low-energy 
fermionic action that rests on the SU(2)-invariant charac- 
ter of the Hamiltonian and a consistent scheme of coher- 
ent states within a functional-integral formulation. We 
show that U(l) and SU(2) gauge fields (the collective high 
energy modes in the SC system) take over the task which 
was carried out by phonons in BCS superconductors and 
play the role of the "glue" that is responsible for the 
formation of the electron pairs. In this sense the present 
work charts a route from the microscopic Hubbard model 
on the square lattice to an effective lower energy action 
that exhibits pairing potential. The paper is organized as 
follows. Section II introduces the model and rotational 
invariant formulation. Section III describes charge and 
spin gauge transformations of fcrmions, which results in 
the phase-angular representation of strongly correlated 
electrons . Section IV is devoted to the evaluation of 
pairing interaction, while Sec. V discusses the effective 
fermionic action. We conclude with Sec. VI. Appendixes 
A, B and C contain miscellaneous results that pertain to 
the technical aspects of the work. 



II. HUBBARD MODEL IN THE ROTATING 
REFERENCE FRAME 

The basic physics of strongly correlated electronic sys- 
tems is the competition between the two tendencies of 
the electron to spread out as a wave and to localize as a 
particle, combined with magnetism. That is, the inter- 
play of the spin and the charge degree of freedom is the 
central issue. These features are encoded in by the Hub- 
bard Hamiltonian - the simplest yet nontrivial model for 



strongly correlated electrons. The relevance of this model 
for superconducting cuprates originates from the obser- 
vation that the one-band Hubbard model tries to mimic 
the presence of the charge-transfer gap of cuprates by 
means of an effective value of the Coulomb repulsion. 
Thus, our starting point is the purely fermionic Hubbard 
Hamiltonian in the second-quantized form 

n = -t J2 [4(r)co(r')+H.c.]+^C/nt(rH(r).(l) 

(rr'),a r 

Here, (r, r') runs over the nearest-neighbor (n.n.) sites, t 
is the hopping amplitude, U stands for the Coulomb re- 
pulsion, while the operator cj, (r) creates an electron with 
spin a =t (= 1); i (= 2) at the square lattice site r. Fur- 
thermore, n(r) — 7i-f(r) +ni{r) is the number operator, 
where na{r) — cJj(r)cQ(r). Usually, working in the grand 
canonical ensemble a term —fi n(r) is added to Ti. in 
Eq.((T]) with fj, being the chemical potential . We treat 
the problem of interacting fermions at finite temperature 
in the standard path-integral formalism^*^ using Grass- 
mann variables for Fermi fields, Cairr) depending on the 
"imaginary time" < t < /3 = l/ksT (with T being 
the temperature) that satisfy the anti-periodic condition 
Ca(rT) — —Ca{rT + /3), to write the path integral for the 
statistical sum Z ^ J [VcDc] e^'^I'^^'^l with the fermionic 
action 

S[c, c] = Sb[c, c] + / drnic, c], (2) 
Jo 

that contains the fermionic Berry term 

Sb[c,c]=^ dTCa{rT)drCa{rT). (3) 

For the problem under study it is crucial to construct a 
covariant formulation of the theory, which naturally pre- 
serves the spin-rotational symmetry present in the Hub- 
bard Hamiltonian. For this purpose, the density-density 
product in Eq.(IT]) we write, following Ref. 21, in a spin- 
rotational invariant way 

= {i^'M - [^^M • S(rr)]'| , (4) 

where S"'(rr)^ = hT,aa' ci(rT)a^^,Ca'irT) denotes 
the vector spin operator (a = x,y,z) with it" be- 
ing the Pauli matrices. The unit vector r2(rT) = 
[sinz?(rr) cos (p{vt), sin ??(rr) sin iy9(rT), cos ??(rr)] written 
in terms of polar angles labels varying in space-time 
spin-quantization axis. Thus, the Hubbard Hamiltonian 
should not change its form under a rotation of the spin- 
quantization axis. This is not apparent in the standard 
form of the interaction in Eq.lU]). The spin-rotation in- 
variance is made explicit by performing the angular inte- 
gration over r2(rT) at each site and time. By decoupling 
spin- and charge-density terms in Eq.(I3|) using auxiliary 
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fields £»(rT) and iV{rT), respectively, we write down the 
partition function in the form, 



Z = y" [Vri] J [VVVq] j [VcDc] 



with Un = 27rn//3 (n = 0, ±1, ±2) being the (Bose) Mat- 
subara frequencies. Now, we introduce the U(l) phase 
field 0(rT) via the Faraday-type relation. 



^-S[n,V,g,c.c] 



(5) 



d(l>{rT) ~ 
(rr) = = V{vt). 



(8) 



where[2?f2] = ]^^^^ £Hi£i£Iii££LpMItl is the spin- 
angular integration measure. The effective action reads 
as 



r0 

5[n,y,e,c,c] = ^ / 



U 



U 



iV{rT)n{rT) + 2g(rT)r2(rr) • S(rT)] 



Since the honiotopy group 7ri[U(l)] forms a set of inte- 
gers, discrete configurations of (/"(rr) matter, for which 
(/)(r/3) - (/)(rO) = 27rm(r), where m(r) = 0,±1,±2,... 
Thus the decomposition of the charge field y(rr) con- 
forms with the basic m = topological sector since 
JI^ 0(rT) = V{rT) = 0. Furthermore, by perform- 
ing the local gauge transformation to the new fermionic 
variables /^(rr). 



drHt [c, c] . 


(6) 










_ Cq (rr) 





z(rT) 
z(rT) 



Urr) 



(9) 



We would like to stress that the fermionic fields in Eq. 

are the physical ones and not due to an enlargement 
of the Hilbert space like in a slave-boson treatment of the 
t-J modelj-i^ As we see in Secs.IIIA-IIIC, the gauge fields 
will arise here by relating an SU(2) rotation in spin space 
and a vector on the two sphere (52). 



III. CHARGE AND SPIN GAUGE 
TRANSFORMATIONS OF FERMIONS 

The interaction term of the Hubbard Hamiltonian is 
quartic in fermion operators. This is a nonlinear problem 
which is not solvable except in some very special cases, 
such as one-dimensional systems. Thus, a standard ap- 
proach is the mean-field approximation, sometimes called 
Hartree-Fock (HF) approximation, in which the quartic 
term is factorized in terms of a fermion bilinear times an 
auxiliary field, which is usually treated classically. How- 
ever, HF theory will not work for a Hubbard model in 
which U is the largest energy in the problem. One has to 
isolate strongly fluctuating modes generated by the Hub- 
bard term according to the charge-U(l) and spin-SU(2) 
symmetries. 



A. U(l) charge transformation 

We swich mow from the particle-number representa- 
tion to the conjugate phase representation of the elec- 
tronic degrees of freedom. To this aim the second- 
quantized Hamiltonian of the model is translated to 
the phase representation with the help of the topolog- 
ically constrained path-integral formalism. To this end 
we write the fluctuating "imaginary chemical potential" 
iV{rT) as a sum of a static Vb(r) and periodic function 
V{rT) = Vb(r) -I- F(rr) using Fourier series, 



V{rT) 



n— 1 



-t- c.c. 



(7) 



where the unimodular parameter |0(rr)p — 1 satis- 
fies z(rr) = e"^^'"^\ we remove the imaginary term 



i dry(rr)n(rr) for all the Fourier modes of the V{yt) 
field, except for the zero frequency. 



B. SU(2) spin transformation 

In the above description, we focused on the charge de- 
gree of freedom of the electron. However, the electron 
has one more degree of freedom being the spin. The 
spin dominates the magnetic properties. The subsequent 
SU(2) transformation from fai^r) to ha{vT) variables, 



(10) 



with the constraint |Ci(rT)p + |C2(rT)p = 1 takes away 
the rotational dependence on r2(rr) in the spin sector. 
This is done by means of the Hopf map. 



" /i(rT) ■ 




■Ci(rr) 


-C2(rr) ■ 




■/ii(rr) 






.C2(rr) 


Ci(rT) _ 




h2{vT) 



R(rr)<T^Rt(rr) =<3--f2(rT) 



where 



R(rr) 



Ci(rr) -C2(rr) 
C2(rT) Ci(rr) 



(11) 



(12) 



that is based on the enlargement from two sphere 5*2 to 
the three sphere 5*3 ^ SU(2). The unimodular constraint 
can be resolved by using the parametrization 



Ci(rr) 



g-i/2[v(rT)+x(rT)] j,Qg 



gi/2[i^(rT)-x(rT)] gjj^ 



2 

t?(rr) 



(13) 



with the Euler angular variables <y5(rr), d{YT) and x(r'''), 
respectively. Here, the extra variable x(i"''") represents 
the U(l) gauge freedom of the theory as a consequence 
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of 5*2 — !■ S3 mapping. One can summarize Eqs. ([9]) and 
(|10p by the single joint gauge transformation exhibiting 
electron operator factorization 

Cairr) = ^ z(rT)i?aa'(rr)/ia'(rr), (14) 

where R(rr) = e-'^-virr)/2^~iay^irT)/2^~ia,x(rr)/2 ^ 
unitary matrix which rotates the spin-quantization axis 
at site r and time r. Equation p^ reflects the com- 
posite nature of the interacting electron formed from 
bosonic spinorial and charge degrees of freedom given 
by i?QQ'(rT) and -z(rT), respectively, as well as remaining 
fermionic core part ha{rT). 

C. Fermionic sector 

Anticipating that spatial and temporal fluctuations of 
the fields Vo(r) and g{rT) will be energetically penal- 
ized, since they are gaped and decouple from the angular 
and phase variables. Therefore, in order to make further 
progress we subject the functional in Eq.® to a saddle- 
point analysis with respect to non-fluctuating (static) 
fields and variables that fluctuations cost energy of the 
order of U. The expectation value of the static (zero- 
frequency) part of the fluctuating potential Vo{r) in the 
charge sector we calculate by the saddle-point method to 
gives 

Voir) = i (^H - ^Uh^ = ip, (15) 

where /2 is the chemical potential with a Hartree shift 
originating from the saddle-point value of the static 
variable Vo{r) with Uh — Uh] -f rihi and Uha = 
(ha{vT)ha{vT)) . Similarly in the magnetic sector, a 
saddle-point evaluation of p(r) reproduces the conven- 
tional HartreeFock equations for a commensurate anti- 
ferromagnet 

p(rr) = (-l)--A, (16) 

where Ac = U{S^'{yt)) sets the magnitude for the Mott- 
charge gap Ac ~ U/2 for U/t ^ 1. The staggeriza- 
tion factor (—1)"" breaks the translation invariance by 
one site which remains by two sites. The term is readily 
handled by going to the reduced Brillouin zone^ Note 
that the notion antifferomagnetic here does not mean an 
actual long-range ordering - for this the angular spin- 
quantization variables governed by the rotational sym- 
metry have to be ordered as well. To summarize, the 
fermionic sector is governed by the effective Hamiltonian 

Hj^,0-E^-(-l)""[^TM/iTM-/ii(rT);ixM] 

r 

-t z(rr)z(r'r)[Rt(rT)R(r'r)]^^/i„(rT)/i,(r'T) 

(rr') ,Q7 

-liYha{rT)ha{vT), (17) 



The chief merit of the gauge transformation in Eq. (fM]) is 
that we have managed to cast the Hubbard problem into 
a system of h fermions submerged in the bath of strongly 
fluctuating U(l) and SU(2) gauge potentials (minimally 
coupled to fermions via hopping term) which mediate the 
interactions. 



IV. PAIRING INTERACTION 



It is well known that the crucial point of BCS theory 
is the existence of an attractive interaction among elec- 
trons, where that phonons play the role of the "glue" that 
is responsible for the formation of Cooper pairs. Here, 
by integration out the bosonic scalar filed that represents 
phonons in the fermionic Hamiltonian, an effective at- 
tractive potential emerges. Now we show that U(l) and 
SU(2) emergent gauge fields (the collective high-energy 
modes in the Hubbard system) take over the task which 
was carried out by phonons in BCS superconductors. In 
a way similar to phonons these gauge fields couple to 
the fermion density-type term via the amplitude i, see 
Eq.Cll, 



~t z(rr)z(r'T)[Rt(rT)R(r'r)]^^/ic(rr)/i^(r'r) 

(rr') ,a7 

(18) 

Thus, in order to obtain an effective interaction among 
fermions we have to integrate out all the bosonic modes 
given by z{tt)^ z(r'r) and Rt(rT), R(r'T) fields. A major 
difference with respect to the BCS theory is that the 
variables to be integrated out are of tensorial nature since 
SU(2) modes carry spin index. To explicitly evaluate 
the effective interaction between fermions by tracing out 
the gauge degrees of freedom, we resort to the cumulant 
expansion To this end we write the partition function 
as Z = J\Dh'Dh]e~^^^'^\ where the effective fermionic 
action is 



Seff [h, h] = -ln [V(j)Vn] e-s[^^^M (19) 



The expression Eq. ([19]) generates a cumulant series when 
expanded with respect to the hopping variable t. The 
relevant second-order term that contains the quartic 
fermionic term becomes 
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X EE ( [R^(rir)R(rlr)]_, [Rt(r2r')R(r^r')] ^ ^ h^{r^r)h^,{r[T)h,ir,r')h,,{r',r') (20) 

/ / SU(2) 



where 



'U(l) 



im 



is the averaging over U(l) phase field while 



-S[0] 



(21) 



(22) 



is the averaging over spin-angular variables. To proceed 
with the evaluation of the effective fcrmion-fermion in- 
teraction one has to develop procedures for effectuating 
both averages which involves calculation of the effective 
actions 5 [(/)] and S [D,] in charge and spin sectors, respec- 
tively. 



A. U(l) average 

The averaging in the charge sector is performed with 
the use of the U(l) phase action (see Appendix A). 



r. Jo 



U 



(rr) 



(23) 



that contains both the kinetic and Berry terms of the 
U(l) phase field in the charge sector. For the U(l) aver- 
age in Eg. ((20)) we get 



z{riT)z{r[T)z{r2T')z{r'2T') j^^^^ W {Sr^r[Sr2r'^ 

+^r,r,V,.)e-[l^-^'|-(-')^H. (24) 
Specializing to the low-temperature limit 



lim 



2 2e'fU 
2 



-^,(25) 



we obtain the result for the U(l) phase average. 



B. SU(2) average 

1. Spin-angular action 

The calculation of the SU(2) average is done with help 
of the effective action that involves the spin-directional 



degrees of freedom O, which important fiuctuations cor- 
respond to rotations. This can be done by integrating 
out fermions 2 = J [DCl] e~'^[^l where 



= - In / [DcpVhVhje 



-S[ip,cl>,^,h,h] 



(26) 



generates the low-energy action in the form S[d] = 
Sq [CI] + Sb [^] + Sj [CI] . The interaction term with the 
spin stiffness becomes 



Sj[n] = ^^-^ V f dTCt{vT) ■ fl{r'T) 



(27) 



with the entiferromagnetic (AF) exchange coefficient 

J(A,) = — (n^-nx)^^ — f-^j . (28) 



it is evident that for [/ — s- oo one has 



From Eq. 

J(Ac) ^ ^ since ^ ^ 1 in this hmit. Thus, m 
the strong-coupling limit, the half-filled Hubbard model 
maps onto the quantum Hcisenberg model. In this limit 
the fermions are bound into localized spins. There is 
no motion of fermions, since they are suppressed by the 
gap for charge fiuctuations. In general the AF-exchange 
parameter persists as long as the charge gap Ac exists. 
However, J(Ac) diminishes rapidly in the U/t ^ weak- 
coupling limit. Because the gauge field is the phase fac- 
tor arising in the inner product of quantum-mechanical 
states - the so-called connection in mathematical lan- 
guage - it is intimately related to the Berry phase term 
Sb [Ci] in the effective action. If we work in Dirac "north 
pole" gauge xi^''') — ~'P{^''') oiis recovers the familiar 
form 



r'^ 

= - V / dTLp{vT)[l - cosi?(rT)]. 



(29) 



Here, the integral on the right-hand side of Eq. ([29|1 
has a simple geometrical interpretation as it is equal to 
a solid angle swept by a unit vector Cl{{), ip) during its 
motion.^ The extra phase factor coming from the Berry 
phase requires some little extra care, since it will induce 
quantum-mechanical phase interference between config- 
urations. In regard to the nonperturbative effects, we re- 
alized the presence of an additional parameter with the 
topological angle or so-called theta term 



A, 
U 



(30) 
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that is related to the Mott gap. In the large-?/ limit, 
one has — > U /2, so that ^ ^ ^ relevant for the half- 
integer spin. The kinetic- energy term in the spin sector 
becomes 







1 



/o 4£ 
-I- 2(p(rT)x(rr) cosi?(rT)] } , 



where 8s — ^/ (2xt) and 



1 



XT 



8J 
I I ft 



(31) 



(32) 



is the transverse spin suusceptibility. 



2. CP^ representation 

It the CP^ representation, the spin-quantization axis 
can be conveniently written as 



(33) 



As a consequence, all the terms in the spin action can be 
expressed as functions of unimodular (rr) variables 
instead of angular variables, which are more complicated 
to be handled. Finally, the action assumes the form 

S [C, C] = E '^^ I 2xtC (rr) • C (rr) 
- {^ly [C (rr) • C (rr) - C (rr) • C (rr)] } 



J} / dTAirTr'T)A rrr r 

/„_,v Jo 



(34) 



(rr') 

with the bond operators 

A (rrr'r) A (rrr'r) ^ (rr) • fl (r'r) + ^ 

. / / ^ Ct (rr) Ci (r'r) - Ci (rr)Ct (r'r) , , 

V2 



5. Canonical transformation of CP^ variables 

In order to achieve a consistent representation of the 
underlying antiferromagnetic structure, it is unavoidable 
to explicitly split the degrees of freedom according to 
their location on sublattice A or B. Since the lattice is 
bipartite, allowing one to make the unitary transforma- 
tion 



CT(rr) 
Ci(rr) 



-Ci(rr) 
CT(rr) 



(36) 



for sites on one sublattice, so that the antiferromagnetic 
bond operator becomes 

^(rrr'r) ^ ^' (rrr'r) = ^ (37) 

a— 1 ^ 

This canonical transformation preserves the unimodular 
constraint of the CP^ fields. Biquadratic (four-variable) 
terms in the Lagrangian cannot be readily integrated in 
the path integral. Introducing a complex variable for 
each bond that depends on "imaginary time"(5(rrr'r), 
we decouple the four- variable terms ^'(rrr'r)y^' (rrr'r) 
using the formula 



- E /^^(^IQI'+QC-C+OC-C) 



I?2Q ^ Y{ d^QivTv'T), 

(rr')r 



(38) 



where (PQ = dReQdhnQ. In a similar manner by intro- 
ducing a local real field A (rr) , we can decouple the sec- 
ond term in the right-hand side, in Eq. p4p . To handle 
the unimodularity condition, one introduces a Lagrange 
multiplier A^(r). Then with the help of the Dirac-delta 
functional, 



Eic(^ 



N 



VXr 



2m 



(39) 

where the variables CT(rr) and Ci(rr) are now uncon- 
strained bosonic fields. Thus, the local constraints are 
reintroduced into the theory through the dynamical fluc- 
tuations of the auxiliary field, so that the statistical 
sum becomes 



X e 



(40) 



where 



+ QC-C + QC-C]. (41) 

Furthermore, by evaluating saddle-point values of the Q, 
a, and fields 

Q.p (rrr'r) = _^(C(rr) • C(r'r)), 



1 



(C(rr) . C(rr)) 



(42) 

and by assuming the uniform solutions Qsp (rrr'r) = Q, 
flsp (rrr'r) = a, and Ai^sp (r) = A^, we obtain for the 
Hamiltonian in the spin-bosonic sector 

W [C, C] = ^ E ^C- (1^^") ^COk i^n) Ac. (ka;„) 



kncT 



(43) 
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with 



where: 



_ Co- (k,cj„) 

_Co- (k - TV,Un) 
C (-k + TT, -W„) 



(44) 



and 



2i0Lu,, -2Q^k r' + ^C 

(45) 

Decouphng of the bond operators in the kinetic term of 
the spin action in Eq. p4|) leads to additional field Q, 
which value is determined from the equation, 



J(A) 

zN 



(46) 



and 



= Ac ± 2gek 



(49) 



(50) 



with = ■^[cos{kx) + cos(fcj^)] as the two-dimensional 
lattice structure factor. 



C. Fermionic action 



while the constraint parameter is the solution of the 
equation 



with 



1 



■/(A) 

TV 



(47) 



Wk 



coth [(3E;^ )] + coth [PE+^ )] 



coth [pe;^ + coth [/3i?+ (^+) 



(48) 



4\/6i2 + > 



We start the calculation of the effective fermionic ac- 
tion with the first-order term with respect to the hopping 
element t 

5«=-i^(z-(rr)z(r'r))u(i) 

(rr') 

[i?t(rr)i?(r'r)]^^)^^^^^ h^{rT)hp{r'T) (51) 

The evaluation of the average with rotational matrices 
gives 



Q/3 



SU(2) 



.5](C4rr)C4r'T))gy(2) h{^T)h^i^'T) = ^ (Co(rT)Co(rV)) [h^{vT)h^{v'r) ~ h^{vT)h^{v'T)] (52) 

a a 

I 



The first-order action is then in the form 



•5f^=-^~E f[h{^r)h^{^'T)-h^{vT)h,{v'T)l (53) 



(rr') 



where t — tgc{A)gg{A), is the renormalized 
hopping, with (7c(d) = (z(rT)z(r'T))u^^^ and 
5s(d) = (CQ(i"r)Ca(r'T))gu(2) being the Gutzwiller- 
type charge and spin renormalization factors. 



D. Second-order contribution to the fermionic 
action 



The calculation of the second-order contribution to 
the effective fermionic action in Eq. ipO]) is more involved 
since the SU(2) averages contain tensorial quantities of 
the form 



M 



OlOl' ,77' 



(rr, rV|rV, rr) 



[Rt(rir)R(r;r)]^^, [Rt(r2r')R(rir')],,,) J54) 
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The sublattice transformation of the CP^ variables in 
Eg. ([55)1 translates to the transformation of the rotation 
matrix R(rr) — » R(rT) matrix 



R(rT) = {iay)R{r'T), 



(55) 



where R(rr) is the transformed form of the rotation ma- 
trix 



R(rr) = 



-C2(rT) -Ci(rT-) 
Ci(rr) -C2(rr) 



(56) 



It is convenient to define the following bond operator 
constructed from the CP^ fields: 

^(rrrV) = ^"^(^^^^^ (^'-^ . (57) 

v2 



With the dedinition in Eq. (l57|) the matrix 



-l^aa' ,77' 

form as 



(rr, r'rjr'T, rr) will be written in a compact 



Maa',^y{rT,r'T\r'T,rT) = 
TT_ TA -TA TT 

ff tA -TA tt 



eta', 77' SU(2) 



where aa' = {11,12,21,22}. Now, we can rewrite 
the second-order fermionic action taking into account the 
non- vanishing averages over CP^ fields (see Appendix C) 
to get 



= _ / dr V Afii,ii(rT,r'r|r'r,rT) V/i„(rT)/i„(r'T)/i„(r'T)/i,(rr) 

" -^0 (rr'> 



+Mii,22(rT, r'r|r'r, rr) 
+Mi2,2i(rT, r'rlr'r, rr) 
+Mi2,i2(rT, r'rlr'r, rr) 



^ ha{rT)ha{r'T)hp{r'T)hp{rT) - ^ ha{YT)ha{v' T)ha{v' T)ha{vT) 



af} 



^ ha{YT)hp{Y' T)hi3{Y' T)ha{YT) - ^ ha{YT)ha{Y' T)ha{Y' T)ha{YT) 



a(3 



^ ha{YT)hp{Y' T)ha{YT')hi3{YT) - ^ ha{YT)ha{Y' T)ha{Y' T)ha{YT) 

a/3 a 



(59) 



In deriving the above result, we made the observation 
that the dynamics of spin variables is slower as compared 
to the charge counterparts, allowing to treat SU(2) vari- 
ables as local intimeR(rr') = R(rr) + (r' — r)9^R(rr) + 



0[(r' — r)^] which leads to nonretarded interactions.— 
Furthermore, with the help of the operator identities 
from Appendix C we can reduce Ea. (|59p to a compact 
form 



S^'-^h I dTj2 [lin{YT)niY'T) + 72A(rrr'r)^;,(rrr'r) + 73S/,(rr) • S,,(r'r) -t- 74n(rr)] , (60) 



where the interaction coefficients 

71 = f{0) + 25^(0) + g'id) + 4/2(d) > 0, 

72 = -2[6/2(d) + 2/2(0)] <0, 
73-4[/2(0)-5'(d)], 

74 = 2g2(d) + 2/2(d)-t- 4^2(0) >0, (61) 
are given in terms of the CP^ normal (g) and anomalous 



(/) correlation functions 

g{Y - r') = - (Ca(rr)Ca(r'r))gy^2^ , 
/(r-r') = (Ca(rr)Ca(r'r))su(2) (62) 

which can be readily evaluated using the propagator of 
the C fields in Eq. (|il)l . 
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FIG. 1: (Color online) Pairing interaction 72 normal- 
ized to the hopping parameter t (upper curve) and the 
antiferromagnetic-exchange parameter J — At'^/U (lower 
curve) as a function of the Coulomb interaction U/t calcu- 
lated at zero temperature and half filling fi = for the two- 
dimensional Hubbard model with nearest-neighbors hopping. 



HAMILTONIAN WITH PAIRING TERM 
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From the result in Eq. (|59p we can deduce the spin- 
singlet pairing possibility in the fermionicsector. To bring 
the kinetic-energy term to a standard form, one performs 
a rotation of the fermionic variables on one of the sublat- 
tices in a manner similar to the bosonic transformation 
in Eq. 



(63) 



As a result the hopping term assumes the conventional 
form that is diagonal in the spin indices 



- 

S^''>[h,h] = -i J2 / K{rT)h^ir'T), 
I /\ "'0 

(rr ) ,a 

while the second-order term is given by 

S(^)[h,h] = y2 dr [7in(rT)n(r'T) 

(rr'> 

-72A(rTr'r)A(rrr'T)] , 



(64) 



(65) 



where 



A(rrr'T) 
A(rrr'T) 



h^{vT)hi{v'T)-hiivT)h^{v'T) 

h{v'T)h^{vT)-h^{v'T)h^ivT) 

V2 



are the bond operators relevant for a singlet pairing. The 
rotational invariance of the right-hand side in Eq. (|65p is 
manifest since 

-Ah{rTY'T)Ah{rTY'T) = 

Sh{vT) ■ S„(r'T) - -^nh{vT)nh{v'T). (67) 

The coefficients 71 and 72 are given by Eq. (pTjl . By not- 
ing that g{d) = and /(O) = one obtains 



71 



4^2 r 



f{d) + -g\0) 



72 = ^ [3/'(d)] - J 



J2(A) 



(68) 



(66) 



where the /, ^-correlation functions can be computed 
with the help of the CP^ propagators: see Eq. (|15)) . The 
effective nonretarded interaction containing 72 in front 
of the A{yty' t)A{vtv' t) term is negative and therefore 
constitutes the attractive potential for fermion pairing. 
We can see that the coefficient 72 is not just given by the 
bare AF exchange J = At'^/U but is renormalized down- 
ward by the quantity Q that is related to the antiffero- 
magnetic spin stiffness as delineated in Sec.IVD dealing 
with the SU(2) spin sector. We have calculated Q self- 
consistently using Eq. fi^ . The result is plotted in Fig.l. 
Note that the pairing interaction survives in rather nar- 
row range of the Coulomb interaction 1.17 < U/t < 1.41. 
This result suggests that superconductivity in the Hub- 
bard model, if possible, represents a rather delicate bal- 
ance between kinetic energy and Coulomb interaction. In 
this context, we note that the deleterious effects of the 
Coulomb interaction superconductivity in cupratcs have 
been largely ignored in the literature. Furthermore, since 
71 > in the density-density term in Eq. (|65[) , many sorts 
of the charge-ordered states can be stabilized, including 
e.g., charge-density wave states, which in general com- 
pete with superconductivity. This is in contrast to the 
BCS theory where the only instability of a Fermi liquid 
is the Cooper instability: the superconducting order is 
generic. 



VI. CONCLUSIONS AND PERSPECTIVES 

The basic physics of strongly correlated electronic sys- 
tems is the competition between the two tendencies of 
the electron to spread out as a wave and to localize as 
a particle combined with magnetism. That is, the in- 
terplay of the spin and the charge degree of freedom 
is the central issue. While there is a growing consen- 
sus that superconductivity in the cuprates arises from 
strong short-range Coulomb interactions between elec- 
trons rather than the traditional electron-phonon inter- 
action, the precise nature of the pairing interaction re- 
mains controversial. While the principal focus of the 
present work is theoretical, the choice of model and the 
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approach is motivated in the experimentally observed 
properties of cuprates. Therefore, in the present work we 
hope to shed some light on this controversial issue with 
the purpose to understand better the physical properties 
of most common model for cuprates. To this end, we 
have discussed in the present work the Hubbard model 
in the spin-rotational invariant formulation which ob- 
serves the important symmetries involved. We presented 
a field-theoretic description of a microscopic model that 
reveals an intimate relationship between the spin-SU(2) 
and charge-U(l) symmetry and pairings. We found that 
the maximal strength of the effective pairing interaction 
parameter is observed in a rather narrow range of U /t 
with the kinetic energy comparable to the Coulomb in- 
teraction. Moreover, the form of the effective fermionic 
action suggests that other competing ordered phases can 
occur simultaneously, which can quench the supercon- 
ductivity substantially. Therefore, the issue of pairing 
interaction is not settling the question about the long- 
range superconducting order in the Hubbard model. As 
far as modeling of cuprates is concerned there is also a 
problem of interplane interaction, entirely omited in the 
present work, which can affect the bulk superconductivity 
considerably. In closing we note that the pairing inter- 
action itself cannot be measured directly: one needs to 
analyze key experiments which reveal fingerprints of it. 
Thus, the continuing experimental search for a pairing 
glue in the cuprates is important and will play an essen- 
tial role in determining the origin of the high-Tc pairing 
interaction. 
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for the Bosonic phase variables is defined as 

g,{vTv'T') = {z{vT)z{v'T')). (69) 

The averaging in this definition is over the U(l) phase 
field and 



(...) 



/ D(/)...e 



J D(f)e-Sol4'] ■ 



(70) 



Here the complex variables z(rT) arc defined as z{rT) = 
gi0(rT)_ rpj^g variables 0(rT) satisfy the following bound- 
ary conditions: 



(r/3) - (f){rO) = 27rm(r). 



(71) 



It is very convenient to satisfy the boundary conditions 
by decomposing the phase field in terms of a periodic 
field (f){rT) and a term linear in r. We set 



(rr) = 4>{rT) + ^m(r) 



(72) 



with (j){(3) = <^(0). Summing over the phase field 
means to integrate all 4'{rT) configurations and perform 
the summation over the integers n. Then we write the 
phase variables <p{rT) in the Fourier-transformed form 

4>{rT) = Me) + 1 g [Mr)e"^"^ + (r)e-^-"1 . 

n=l 

(73) 

The weight of the averaging in the expression of the phase 
correlator is given by the following exponential: 
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VIII. APPENDIX A:U(1) PHASE AVERAGES 



-So[4>] = lo^ d.T^^{rT)/U^ -So[m] 



(74) 



where the action 5*0 [m] is the topological part of the ac- 
tion and is given by 



So[m] = 



r n 
-2^2 



I3U 



U 



m(r) 



(75) 



In this section we evaluate the expression for U(l) 
phase propagator. Two point phase-phase propagator 



Now we evaluate the average. We write first the non- 
topological part of the action 



J[D(f)]e-y^u^r^n\'l>n{r)\^ 





' up ■ 

4a; 




uT) sin(a'nT')]^-t-[cos(uJn'r)— COs(uJnT')]^| 




\ Ul3] 





(76) 
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By using the identity 

[sin(x) — sm{y)]^ + [cos{x) — cos{y)]^ 
= 2 — 2[cos(a;) cos{y) + sin(x) sin(?/)] 
= 2 — 2cos(x — y) 



5; 



.(rrrV) = S,,,e-^^<\^-^'^-<^^-^'^'/^] 



(80) 



(77) 



one obtains 

5,(rTr'r') = ^^r'^-^/'^^-i iK{i-cosK(r-r')]}. (78) 

Now, in order to calculate the sum in the exponential we 
use the following identity: 



Now we are ready to calculate the four-point phase cor- 
relator: 

(z(riT)z(r;)z(r2T')z(r2T'))u(i) = 

(z(riT)z(r;))u(i) (z(r2r')z(r2T'))u(i) + 

+ (z(riT)z(rir'))u(i) (z(r;T)z(r2r')>u(i) (81) 



4 4cos(a;„r)x 



(79) 



By using the result in Eq. HO]) we get Eq. (|24l) . 

IX. APPENDIX B:SU(2) AVERAGE 



where —P<x<f3. And finally we get 



J 



The composition formula for the rotational matrices in the angular representation is given by 



1 



i?t(rr)i?(fr) = 



yi + »(rr)ri(rV) exp(|$) ^1 - ri(rr)0(rV ) exp(|^) 
VI - fl(rr)rj(r'r) exp(-f $) ^1 + rj(rr)J7(r'T) exp(-|$) 



(82) 



where $ = $[r2(rr), r2(r'T), z] is the signed solid angle spanned by the vectors Jl(rT), r2(r'r) and z with $ 
$[r2(rr), — r2(r'r)] — 2(p{rT). In the complex projective representation, Eq. ([82|) reads 



i?t(rr)i?(rT) = 



Ci(rT)Ci(r'T) + C2(rr)C2(r'r) -Ci(rT)C2(r'T) + C2(rr)Ci(r'T) 
-C2(rr)Ci(r'T) + Ci(rT)C2(r') C2(rT)C2(r'T) + Ci(rT)Ci(r'T) 

I 



(83) 



The form of Eq. ((83| suggests the use of the bond opera- 
tors defined by Eqs. ((35|l and ([57|l . so that the product of 
rotational matrices can be written in a compact form 



R^fTrT)R{rT) = 



T -A 
A T 

T A 

-A T 



(rrr'r) 

(rrr'r). (84) 



Under the transformation the bond operators become 
-Ci(rT)C2(r'r) + C2(rr)Ci(r'T) 



With the help of the above equation it is easy to write 
down the components of the M matrix 



J 



V2 

A{vtv't) ~* A!{vtv't) 

Ci(rr)C(rV)+C2(rr)C2(rV) 
V2 



(85) 



Afi2;2i(rr,r'r|r'T,rr) = ( [i?t(rr)i?(rV)] {R\v' t)R{vt%^) ^^^^^^ 

= (^'(rrr'T)^'(rrr'r))su(2) - ( [Ci(r'r)Ci(rr) + C2(r'r)C2(rT)] [Ci(rr)Ci(r') + C2(rr)C2(r'T)])gy^^^ 
= /5^Co(r'r)Ca(rr)C^(rr)CMr'r)\ 

\a,/3 / su(2) 

Furthermore, by implementing the Wick theorem to the CP^ averages 

(Ca(r'T)C«(rT)C/3(rT)C/3(r'T))gy^2) = (C«(i"'T')Ca(rT))gu(2) (C/3(i-T)C/3(r'T-))gu(2) 

+ (Ca(r'T)C^(rT))gy^2) (C"MC/3(i"'^))su(2) 

+ (Ca(r'T)C^(r'T))su(2) (Ca(rT)C/3(rT))gy(2^ 

= ^QQ<5,5/3/(-d)/(d) 5a/3<5/3Qff(ti)g(-d) -I- jQ/3(5^jag(0).g(0) 



(86) 



(87) 
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In a similar manner 



Afn,n(rT,r'T|r'T,rT) = 2 [^^(O) - f{d)] , 
Mn,22(rr,r'T|r'T,rr) = 2 [f{0) - g^d)] , 
Mi2;2i(rT, r'T|r'T,rT) = 4/2(d) + 2g(d) + 2^(0), 
Mi2;i2(rT,r'T|r'T,rT) = - [6f{d) + 2f{0)] 



where we have used that /(— d) = /(d) and g{—d) = 
g{d). It is not difficult to see that 

Mii,ii(rT,r'T|r'T,rr) = M22,22(rr, r'r|r'r, rr), 
^ii;22(rT,r'r|r'T,rr) = M22,ii(rT, r'r|r'r, rr), 
A^i2;2i(rT, r'r|r'r, rr) = M2i;i2(rT, r'r|r'r, rr), 
Afi2;i2(rT, r'rlr'r, rr) = M2i;2i(rT, r'rlr'r, rr), (89) 

while all the remaining components of the M matrix van- 
ish. 



one can prove the following useful identities that involve 
four-fermion products that appear in the second-order 
cumulant expansion. For the spin and charge products, 
we have 



Shirr) ■ S„(r'T) 



nhir'r) 



Ah{rTr'T)Ah{rTr'T) _ Th{rTr'T)Th{rTr'T) 
2 2 
nh{rT)nh{r'T) _ n/t(rV) 

2 ~ 2 _ 

+AhivTr'T)Ah{rTr'T) - Th{vTr'T)Th{vTr'T). (92) 
Furthermore, for the prpducts of fcrmionic variables that 
appear in the calculation of the M matrix, one finds 



X. APPENDIX C:USEFUL OPERATOR 
IDENTITIES 

By introducing the fermionic representation of the 1/2- 
spin operators 



S/i(rr) = ^ X! ha{YT)&afsh0{YT) 

and the following bond operators: 

M^irTv'r) - ^T(rT)/iT(r'T) + /'4(rr)/ii(r'r) 

AhivTY'r) 



(90) 





V2 


h]{vT)hi{r'T 


)-h^{vT)h^{v'T) 








)+HvT)h^{v'T) 



'^ha{vT)ha{v' T)hi3{Y' T)h(3{YT) = nh{vT) 

nh{vT)nh{v'T) . . „ . 
^ 2Sh(rr) • Sh(r r), 

^ha{vT)hi3{Y'T)hi3{Y'T)ha{YT) = 2nh(rr) 
-n,i(rr)n,i(r'r), 

^ha{vT)hf3{v' T)ha{v' T)hi3{vT) = n/,(rr) 

a/3 

-2^UrTr'T)^UrTr'T), 



(93) 



x/2 



where, by simple inspection, one can find that the spin- 
rotational symmetry is apparent. 
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Epjojo ^ 10'^. Since the reriormalizatiori of the Coulomb 
pseudopotential is logarithmic, this largo value of the re- 
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the other hand, Ef/wd < Ef/2Ao ~ 5, where Ao is the 
superconducting gap, so the "HTSC" materials are clearly 
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